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ABSTRACT 

This  paper  presents  an  algorithm  for  computing  the  optimal 
target  path  for  two  aircraft  traversing  a  target  area  from 
different  directions.   There  are  constraints  on  the  maneuver- 
ability of  each  aircraft  which  prohibit  it  from  attacking  every 
target.   The  algorithm  chooses  a  subset  of  targets  whose  de- 
struction will  yield  maximum  value  to  the  attacking  force. 
The  basis  of  the  algorithm  is  the  branch  and  bound  method, 
with  upper  bounds  computed  by  dynamic  programming.   Several 
variations  are  considered,  such  as  payload  limit,  an  increased 
number  of  aircraft  from  each  direction,  and  a  three-directional 
attack.   An  example  problem  is  solved  using  the  basic  model. 

A  Fortran  IV  computer  program  is  included.   Computation 
time  versus  problem  characteristics  is  discussed. 
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I.   INTRODUCTION 

Optimal  assignment  of  targets  to  aircraft  on  a  strike 
mission  can  greatly  increase  the  effectiveness  of  an  air  attack 
against  ground  installations.   The  purpose  of  this  paper  is  to 
present  an  algorithm  which  will,  given  the  location  and  rela- 
tive military/industrial  worth  of  key  enemy  positions,  select 
a  sequenced  subset  of  targets  whose  destruction  will  yield 
maximum  value  to  the  attacking  force.   It  assumes  that  the 
decision  maker  has  complete  knowledge  of  the  targets  and  their 
value.   A  stochastic  extension  can  easily  be  incorporated 
when  considering  hardened  targets.   This  is  done  by  multiply- 
ing the  probability  of  killing  a  particular  target  by  its  value. 

The  initial  course  of  an  aircraft  approaching  a  target 
area  is  denoted  by  <j>.   Due  to  anti-air  defenses,  there  is  a 
limit  placed  on  the  maximum  number  of  degrees  which  an  aircraft 
may  deviate  from  its  initial  course.   This  angular  deviation  is 
denoted  by  6.   its  effect  is  to  restrict  the  aircraft's  movement 
at  any  point  to  a  cone  whose  vertex  angle,  29,  is  bisected  by 
the  aircraft's  initial  course,  <j>.   Thus  at  any  point  (x,y)  in 
the  target  area,  the  aircraft's  course  options  are  between 
<j>-6  and  <}>+9.   This  creates  an  ordering  among  the  targets  which 
dictates  the  sequence  in  which  they  must  be  considered.   The 
allowable  deviation  is  pictured  in  Figure  1. 


Not  Possible 


Not  Possible 


(x,y) 


Figure  1.   Cone  of  Allowable  Course  Deviations 


Because  of  the  course  constraint,   there  are  only  a  limited 
number  of  paths  through  the  target  area  which  are  feasible. 
The  task  of  the  decision  maker  is  to  choose  that  set  of  tar- 
gets which  yields  the  optimal  combined  value. 

References  [1]  and  [2]  present  a  method  for  determining 
the  optimal  set  of  targets  for  one  aircraft  traversing  an  area 
containing  MN  targets.   The  problem  is  formulated  as  an  MN+1 
stage  dynamic  programming  problem  [Ref .  3] .   The  targets  are 
numbered  in  decreasing  order  from  an  imaginary  line  drawn  tan- 
gent to  the  boundary  at  which  the  aircraft  enters  the  target  area. 
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The  stages  are  represented  by  straight  lines  drawn  through 
the  targets,  parallel  to  the  boundary  tangent  line.   Stage  n 
corresponds  to  the  parallel  line  drawn  through  target  n. 
The  starting  point  of  the  aircraft,  outside  the  target  area, 
is  stage  MN+1.   It  is  located  so  that  any  target  is  accessible 
from  stage  MN+1.   The  stage  diagram  is  depicted  in  Figure  2. 
Should  two  targets  be  equidistant  from  the  boundary  tangent 
line,  either  one  may  be  assigned  the  next  sequential  number, 
and  the  problem  reduces  to  MN  stages. 
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Figure  2.   The  Stage  Diagram  for  One  Aircraft 
Entering  from  below  the  Target  Area. 

8 


It  is  important  to  realize  that  were  the  aircraft  to  enter 
the  target  area,  say,  from  the  right,  the  stage  diagram  would 
appear  as  in  Figure  3 . 
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Figure  3.   The  Stage  Diagram  for  One  Aircraft  Entering  Target 
Area  from  the  Right. 


The  state  variable  x  denotes  the  lateral  position  of  the 

n  c 

aircraft  at  stage  n.   D   is  the  decision  as  to  the  heading  of 
the  aircraft  as  it  moves  from  stage  n  to  stage  n-1.   D   is 
restricted  to  lie  in  the  set  of  feasible  headings  S  ,  from 
<j>-8  to  <j>+9.   The  state  variable  x  is  then  a  function  of  x  ., 
and  D   , .   This  function  is  referred  to  as  the  stage  trans- 
formation t. 

The  return  function  for  stage  n  is  denoted  by  r  .   Letting 
p  be  the  lateral  position  of  target  n,  and  V  be  its  value, 


r  (x  )  =    V      if  x  =  p 
n  n        n         n   ^n 


otherwise 


The  problem  is  then  written  as 

MN 
Maximize    Z    r   (x  ) 
n-1   n   n 


Subject  to:  ^n=t  (xn+]_'Dn+]_)   n=1' ' 


MN 


D„  c  S„  n=2,...,MN+l 

n    n 


An  efficient  algorithm  for  solving  this  problem  is  dis- 
cussed in  Section  VI,  but  the  problem  can  also  be  easily  solved 
graphically.   To  do  so,  the  course  deviation  angle  must  be 
viewed  from  the  perspective  of  the  stage  diagram.   Figure  4 

shows  whether  target  n  is  feasible  for  various  values  of  x  . , . 

n+i 
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STAGE       X"+1 /   X"+1-    \      X"+1 

n+1        not  feasible      feasible       not  feasible 


Figure  4.   Cone  of  Feasibility 

Using  Figure  2  as  an  example,  the  maximum  return  possible 
for  stages  (n-1)  to  1  is  recorded  on  each  stage  line  n, 
n  =  1,  ...,  MN+1,  for  every  possible  value  of  x  .   Assuming 
the  value  of  each  target  is  one,  Figure  5  gives  the  solution 
for  one  aircraft  entering  from  below  the  target  area. 

Tracing  back  from  stage  MN+1,  the  optimal  set  of  targets 
is  8,  7,  6,  4,  3,  2,  for  a  value  of  six. 

Dynamic  programming  can  extend  this  problem  to  M  aircraft 
by  increasing  the  number  of  state  variables  to  M  [Refs.  1,  2, 
and  3] ,  as  long  as  all  M  aircraft  are  attacking  from  the  same 
direction.   But  when  the  aircraft  enter  the  target  area  from 
different  directions,  a  pure  dynamic  programming  solution  is 
no  longer  possible,  since  the  association  between  the  target 
and  the  stage  is  no  longer  valid.   A  target  in  stage  three 
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for  one  aircraft  may  be  in  stage  ten  for  another,  and  the  re- 
cursive equations  have  no  meaning. 


STAGE 

1 

2 

3 

4 
5 


Figure  5.   Graphical  Solution  to  the  One-Aircraft  Problem. 
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This  paper  presents  a  solution  to  the  problem  of  two  air- 
craft attacking  a  target  area  from  different  directions.  A 
model  is  formulated,  and  the  algorithm  developed.  An  example 
problem  is  included.  Also  included  is  a  computer  program  for 
implementing  the  algorithm  and  a  discussion  of  its  effective- 
ness. Because  the  terminology  becomes  quite  involved,  a 
glossary  of  terms  is  provided  in  Appendix  A. 
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II.   PROBLEM  FORMULATION 

In  formulating  this  problem,  it  is  assumed  that  if  an 
aircraft  attacks  a  particular  target,  that  target  is  destroyed. 
In  other  words,  there  are  no  misses.   As  was  stated  in  Section 
I,  this  does  not  preclude  a  stochastic  approach  for  hardened 
targets  where  each  has  a  probability  of  being  destroyed. 

For  ease  of  notation,  the  two  aircraft  are  denoted  air- 
craft A  and  aircraft  B. 

The  targets  are  numbered  as  in  Section  I,  but  now  each 
target  has  two  numbers,  one  with  respect  to  aircraft  A,  and  one 
with  respect  to  aircraft  B.   They  are  denoted  A.  and  B.  ,  respec- 
tively, i  =  1,  ...,  MN,  j  =  1,  ...,  MN.   It  is  critical  to  real- 
ize that  if  any  target  A.  has  the  same  coordinates  as  any  tar- 
get B-,  then  that  A.  is  the  identical  target  B..   In  fact,  for 
every  A.  there  will  be  a  B .  identical  to  it.   Figures  2  and  3, 
which  are  identical  target  areas,  should  clarify  this  point. 

The  value  of  A.  will  be  denoted  by  V[A.] ,  and  the  value  of 
B.  by  V[B.]  . 

A  feasible  set  of  targets  for  a  single  aircraft  is  one  in 
which  the  course  required  to  go  to  each  successive  target  is 
within  the  allowable  course  deviation  cone  described  in  Sec- 
tion I.   PA (I)  denotes  one  of  the  feasible  sets  of  targets 
for  aircraft  A.   It  is   convenient  to  think  of  PA(I)  as  a 
path  of  targets  which  aircraft  A  will  attack.   Any  PA (I)  com- 
pletely disregards  aircraft  B  in  that  it  is  computed  as  if  only 
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one  aircraft  were  attacking.   Similarly,  PB(J)  denotes  one 
of  the  feasible  sets  of  targets  for  aircraft  B. 

The  value  of  a  set  of  targets,  or  a  path,  is  the  summa- 
tion of  the  individual  target  values  comprising  that  set.   The 
values  of  PA(I)  and  PB(J)  are  denoted  by  V[PA(I)]  and 
V[PB(J)],  respectively. 

In  order  for  the  solution  of  a  two-aircraft  problem  to  be 
feasible,  the  sets  of  targets  for  aircraft  A  and  aircraft  B 
must  be  mutually  exclusive,  that  is,  they  must  have  no  targets 
in  common.   This  is  so  because  one  aircraft  is  sufficient  to 
destroy  the  target.   No  additional  return  would  be  realized 
by  the  other  aircraft  attacking  the  same  target.   If  this 
feasibility  constraint  were  not  required  for  optimality,  aircraft 
A  (or  aircraft  B)  might  forego  the  opportunity  to  attack  other 
targets  in  order  to  attack  the  target  that  both  aircraft  have 
in  common. 

The  optimal  solution  is  achieved  when  PA (I)  and  PB(J)  are 
chosen  so  as  to 

Maximize  V[PA(I)]    +V[PB(J)] 

Subject    to:       PA  (I)     C\       PB(J)    =   0. 
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III.   THE  ALGORITHM 

The  branch  and  bound  method  forms  the  basis  of  the  solution 
algorithm  for  the  two-aircraft  problem.   It  is  discussed  in 
theory  in  Refs.  4,  5,  6,  and  7.   Only  a  small  fraction  of  the 
possible  solutions  to  the  problem  is  actually  enumerated.   The 
remaining  solutions  are  eliminated  from  consideration  through 
the  application  of  bounds  that  establish  that  such  solutions 
cannot  be  optimal . 

The  algorithm  begins  by  considering  all  possible  combina- 
tions of  paths  for  both  aircraft.   It  then  breaks  this  set  of 
all  possible  combinations  into  smaller  and  smaller  subsets,  and 
calculates  for  each  an  upper  bound  on  the  value  of  the  best 
paths  contained  therein.   The  bounds  determine  the  partitioning 
of  the  subsets  and  eventually  identify  an  optimal  path  for  both 
aircraft.   The  branch  and  bound  method  represents  the  subsets 
as  nodes  of  a  tree  and  the  partitioning  of  the  subsets  as  a 
branching  of  the  tree. 

Node  one  consists  of  all  possible  combinations  of  paths  for 
both  aircraft.   Using  the  single  aircraft  dynamic  programming 
(D.P.)  method  of  Section  I,  the  optimal  path  for  aircraft  A 
is  computed  for  the  direction  from  which  A  enters  the  target 
area.   This  path  is  denoted  by  PA(1) .   By  the  same  method,  the 
optimal  path  for  aircraft  B  is  computed  for  the  direction  from 
which  B  enters.   This  resultant  path  is  denoted  by  PB(1).   It 
is  important  to  realize  that  both  PA(1)  and  PB(1)  are  computed 
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as  single  aircraft  optimizations,  and  that  all  MN  targets  are 
possible  elements  of  PB(1) ,  even  those  in  PA(1) . 

If  PA(1)  and  PB(1)  have  no  targets  in  common,  then  the 
paths  form  the  optimal  solution  to  the  two-aircraft  problem,  the 
value  of  which  is  V[PA(1)]  +  V[PB(1)],  and  the  algorithm  stops. 
However,  if  there  are  targets  in  common,  then  as  was  mentioned 
in  Section  II,  the  solution  is  not  feasible.   The  summation 
V[PA(1)]  +  V[PB(1)]  represents  instead  an  upper  bound,  denoted 
by  UB(1),  for  the  optimal  solution.   The  summation  of  the 
values  of  the  points  of  intersection  of  PA(1)  and  PB(1)  is  de- 
noted by  VINT(l) .   The  target  in  the  intersection  which  has 
the  highest  value  is  denoted  by  X^  (1)  .   This  target  could  be 
written  in  terms  of  aircraft  A  or  aircraft  B.   When  the  dis** 
tinction  is  necessary,  X   (1)  will  be  written  as  either  X   (1) 
or  X^-CD-.,  for  aircraft  A  and  aircraft  B,  respectively.   Geo- 

graphically,  however,  X^d),  XAB(DA'  and  XAB(1)Bare  identical« 

The  goal  is  to  find  optimal  paths  for  aircraft  A  and  B 
which  have  no  targets  in  common.   The  nature  of  the  D.P.  solution 
for  a  single  aircraft  path  is  such  that  it  seeks  out  those  fea- 
sible targets  with  the  highest  value.   In  node  one,  both  A  and 
B  sought  XAB(1).   If  the  set  of  path  combinations  in  node  one 
was  restricted  so  that  aircraft  A  had  to  take  X-^d),  and  air- 
craft B  could  not  take  X^Cl)-,  the  same  target,  then  that 
point  of  intersection  would  be  eliminated.   But  perhaps  the 
optimal  solution  requires  B,  not  A,  to  take  X-_(l).   It  could 
even  require  that  neither  path  include  XAB(1).   Therefore, 
to  include  all  possibilities,  the  set  of  node  one  is  broken 
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into  two  subsets,  one  which  requires  that  the  path  of  aircraft 

A  exclude  X   (1) ,  giving  aircraft  B  the  option  of  taking  it 

or  not,  depending  on  the  single  aircraft  D. .  P.  solution  for  B. 

The  other  subset  requires  that  the  path  of  aircraft  B  exclude 

X  _(1),  giving  aircraft  A  the  option.   Thus  node  one  branches 
Ad 

to  form  nodes  two  and  three. 

The  restrictions  placed  on  each  aircraft  at  node  I  are  de- 
noted by  R(I) .   If,  for  example,  R(I)  =  A.,  B.,  B,,  the  path 

1    3    Jc 

of  aircraft  A  would  be  required  to  exclude  target  A.,  and  the 
path  of  aircraft  B  would  be  required  to  exclude  targets  B . 
and  B,  . 

The  branching,  with  restrictions,  is  illustrated  in 
Figure  6. 


R(2>=XAB(1)A    /  \   R(3>=XAB(1>B 


Figure  6.   The  Start  of  a  Tree 

Now  consider  node  two  with  the  restriction  vector  R(2) 
containing  the  single  element  xab^A'   A<?a;i-n  using  the 
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single  aircraft  D.P.  method,  only  this  time  with  the  restric- 
tion that  the  path  of  aircraft  A  not  include  target  X^d),, 
the  optimal  path  for  aircraft  A  is  computed.   The  restriction 
can  be  incorporated  into  the  solution  techniques  of  Section  I 
by  temporarily  assigning  a  large  negative  value  to  X   (1)  , 
thereby  making  it  highly  unattractive  as  an  element  of  the 
optimal  path.   Denote  the  resultant  path  PA (2) ,  and  return 
the  original  value  to  XAB(D«   As  was  noted,  restricting 
PA (2)  to  exclude  XaB(l)A  places  no  restrictions  on  aircraft 
B.   Therefore,  PB(2)  will  be  identical  to  PB(1).   Just  as 
with  node  one,  UB(2),  VINT (2),  and  X   (2)  can  be  computed. 
Node  three  is  considered  next,  with  restriction  vector 
R(3)  containing  the  single  element  X__(l)_.   Since  this  means 

Ad    a 

that  the  path  of  aircraft  B  must  exclude  X--,(l)_,  a  large 

Ac     ti 

negative  number  is  assigned  to  X  _(1)  ,  PB(3)  is  computed, 
and  the  original  value  is  returned  to  XAB(1)_.   Since  R(3) 
places  no  restrictions  on  PA (3) ,  it  is  identical  to  PA(1). 
Next,  UB(3),  VTNT(3),  and  X   (3)  are  computed. 

A  terminal  node  is  one  from  which  branching  may  still 
occur.   Nodes  two  and  three  are  terminal  nodes.   Since  only 
two  branches  may  emanate  from  any  node,  node  one  is  no  longer 
terminal,  and  need  not  be  considered  further. 

The  next  step  in  the  algorithm  is  to  choose  the  terminal 
node  which  has  the  highest  upper  bound.   Assume  it  is  node 
J  (at  this  point,  J  is  either  two  or  three,  whichever  has  the 
higher  upper  bound) .   Node  J  then  branches  to  form  nodes  four 
and  five.   The  restriction  vector  R(4)  will  equal  R(J)  with 
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the  addition  of  element  X._ (J) _ .   The  restriction  vector 

AB    A 

R(5)  will  equal  R(J)  with  the  addition  of  element  X__(J)_. 

AB    B 

Again  the  restricted  paths  are  computed  and  the  algorithm 
continues  until  the  stopping  condition  is  met. 

The  branching  is  always  done  in  pairs,  a  left  branch  and 
a  right  branch,  as  was  shown  in  Figure  6.   No  valid  upper 
bound  comparisons  can  be  made  until  both  branchings  have 
been  performed  and  the  upper  bounds  of  both  new  nodes  have 
been  computed. 

The  stopping  condition  for  the  algorithm  occurs  when, 
following  a  double  branching,  a  terminal  node  is  found  whose 
upper  bound  is  greater  than  or  equal  to  all  other  terminal 
node  upper  bounds,  and  whose  paths  for  aircraft  A  and  aircraft 
B  are  mutually  exclusive,  that  is,  they  have  no  targets  in 
common.   This  solution  is  optimal  because  its  value  is  equal 
to  the  upper  bound  of  that  node,  making  it  the  best  solution 
possible  for  that  node.   And  since  its  value  is  at  least  as 
good  as  the  best  solution  possible  for  all  other  terminal 
nodes,  it  is  a  global  optimum. 

The  steps  of  the  algorithm  are  summarized  in  Figure  7 
at  the  end  of  this  section. 

There  is  one  other  calculation  which  can  be  made  at  each 
node  which  is  of  some  interest.   The  lower  bound  for  each 
node  I  is  denoted  by  LB  (I)  and  is  equal  to  UB(I)  minus  VINT  (I)  . 
This  is  a  lower  bound  because  if  all  the  points  of  intersection 
on  PA (I)  and  PB(I)  were  given  to  either  aircraft  A  or  air- 
craft B,  the  value  of  the  resultant  solution  would  be  LB (I) . 
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If  any  terminal  node  I  has  a  lower  bound  on  the  optimal  value 
which  is  greater  than  or  equal  to  the  upper  bound  of  any  other 
terminal  node  J,  then  node  J  may  be  completely  dropped  from 
consideration.   It  is  said  to  be  fathomed,  and  is  no  longer 
terminal.   This  will  not  speed  up  the  algorithm  or  reduce  the 
number  of  branchings  required,  since  the  algorithm  would 
never  branch  from  node  J  anyway.   However,  if  computer  storage 
space  were  critical,  it  would  be  advantageous  to  incorporate 
lower  bounding,  since  once  a  node  was  fathomed  it  could  be 
removed  from  storage.   It  will  not  be  employed  in  this  al- 
gorithm. 

The  algorithm  guarantees  that  an  optimal  solution  will  be 
found.   However,  it  suffers  from  a  limitation  common  to  all 
branch  and  bound  methods.   For  any  untried  problem,  it  is  im- 
possible to  tell  beforehand  exactly  how  much  computation  will 
actually  be  necessary  to  find  the  optimal  solution.   Depending 
on  the  way  the  problem  is  set  up,  it  could  converge  to  the 
optimum  very  quickly,  or  for  a  large,  difficult  problem  it 
could  require  such  excessive  branching  that  it  becomes  compu- 
tationally prohibitive.   This  would  be  the  case  if  the  allow- 
able course  deviation  were  very  large,  the  angle  between  at- 
tackers very  small,  and  there  were  multiple  optimal  paths. 
However,  in  sample  problems  of  target  optimization,  the  al- 
gorithm converges  very  quickly,  as  will  be  shown  in  Section 
VI. 
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the  Algorithm. 


IV.   A  TEST  PROBLEM 

Two  aircraft  are  to  attack  a  target  area.   Aircraft  A 
has  an  initial  heading  of  due  north,  and  aircraft  B  is  heading 
due  west  as  they  approach  the  area.   Each  aircraft  has  an 
allowable  course  deviation  of  forty-five  degrees. 

The  target  positions  and  values  are  given  in  Table  I. 
The  positions  are  given  in  terms  of  the  cartesian  plane,  with 
the  positive  Y  axis  pointing  due  north. 

TABLE  I 
TARGET  POSITIONS  AND  THEIR  VALUES 


Target  Positions  (X,  Y) 

Target  Values 

9,  12 

1 

13,  5 

2 

4,  0 

4 

8,  6 

2 

6,  7 

1 

11,  14 

3 

2,  1 

2 

10,  11 

1 

5,  10 

5 

7,  4 

1 

3,  8 

3 

1,  13 

5 

15,  15 

2 
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Figure  8.   Stage  Diagram  for  Aircraft  A 
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Figure  9.   Stage  Diagram  for  Aircraft  B 
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The  stage  diagrams  for  aircraft  A  and  aircraft  B,  including 
allowable  course  deviation  cones,  are  illustrated  in  Figure  8 
and  Figure  9,  respectively. 

Beginning  with  node  one,R(l)  equals  the  null  set,  meaning 
there  are  no  restrictions  on  the  path  of  either  aircraft. 
Using  the  single  aircraft  D.P.  method  for  aircraft  A,  it  is 
found  that 

PA(-L;  =  A,  .j  ,  A,,,  A~  ,  A, 

PA(1)  =  A, »,  An'  at    Ac  is  also  optimal,  and  either  may  be 
chosen.   For  this  example,  the  former  is  used.   Similarly, 
for  aircraft  B, 

PB(1)  =  ®13 '  ®H'  ^9'  ^5'  ^1 
From  PA(1)  and  PB(1),  the  following  values  are  computed: 

XAB(1)  *  A3  =  Bl 
VINT(l)  =  5 

UB(1)  =  V[PA(1)]  +  V[PB(1)]  =  29 
Since  both  paths  have  target  A3  =  B,  in  common,  they  do  not 
represent  a  feasible  solution  for  the  two-aircraft  problem. 
Node  one  branches  to  form  nodes  two  and  three. 

R(2)  =  XAB(1)A  =  A3 
R<3>  =  XM(1)B  =  Bx 

The  tree  at  this  point  is  illustrated  in  Figure  10. 
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R(2)=A3 


R(3)=B. 


Figure  10.   Start  of  the  Tree 

Considering  node  two,  PA (2)  and  PB(2)  are  computed  sub- 
ject to  the  restriction  that  PA (2)  may  not  include  target  A, 
There  are  no  restrictions  on  PB(2). 

PA  \£.  i     =    A-.  -j  i    A-.  -.  i     A—  /  A  - 
PB(2)  =  B13,  B1X,  Bg/  B5,  B1 

The  target  of  intersection  is  Afi  =  B,.. 

XAB(2)  '   A6  *  B5 
VINT (2)  =  5 

UB(2)  =  29 
Moving  to  node  three,  PA (3)  and  PB(3)  are  computed  sub- 
ject to  R(3)  which  states  that  PB(3)  must  exclude  target  B,. 

FA  \-i  )     =    &■*  ^  ,  A-.  -I  ,  A-» ,  An 
PB(3)  =  B.  ^  i    B,  ,,  ^9'  ^5'  ^3 
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The  target  of  intersection  is  A_  =  B_  . 

XAB(3)  *  A7  =  B3 
VINT (3)  -  3 

UB(3)  =  27 
Since  a  double  branching  has  been  completed,  the  upper 
bounds  of  all  terminal  nodes  must  be  compared.   The  terminal 
nodes  are  nodes  two  and  three.   Node  two  has  the  highest  upper 
bound.   Since  PA (2)  and  PB(2)  have  a  target  in  common,  the 
paths  cannot  be  feasible  for  the  two-aircraft  problem.   There- 
fore, node  two  branches  to  form  nodes  four  and  five.   The 
branching  tree  expands  to  Figure  11. 


R(4)=A3,A6 


R(2)=A 


R(3)=B. 


Figure  11.   Expanded  Branching  Tree 


Node  four  is  subject  to  the  restrictions  from  node  two 
in  addition  to  the  restriction  that  PA (4)  cannot  include  A, 

cA.('i)     —    A^3 ,  A-.  ■■  ,  A- ,  -"c/   4'   O 
PB(4)  -  B13,  BX1,  Bg,  B5,  B1 
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. 


Here  there  are  two  targets  of  intersection,  A.  =  Bg,  and 
A2  =  Bn. 

XAB(4)  =  A2  =  Bn 

VINT  (4)    =   V[A4    =   B9]    +   V[A2   =    B-^]    =    4 

UB(4)  -  28 
Node  five  has  R(5)  equal  to  R(2)  with  the  added  restric- 
tion that  PB(5)  not  include  XAB(2-)«.   Thus  R(5)  =  A, ,  B_. 
PA (5)  must  exclude  target  A,,  and  PB(5)  must  exclude  target 
B5. 

PA  \0)     —    A-.  ~  i     A.  i  ,  A—  /  A r 
PB(5)  =  B13,  B1;L,  Bg,  B1 

There  are  no  targets  of  intersection,  but  this  does  not  mean 
that  the  optimal  solution  has  been  found,  since  terminal  node 
upper  bounds  have  not  yet  been  compared. 

XAB^5)  *  ' 
VINT (5)  ■  0 

UB(S)  =  24 

Having  completed  a  double  branching,  upper  bounds  are 
now  compared  for  terminal  nodes  three,  four,  and  five,  and  it 
is  found  that  node  four  has  the  highest  upper  bound.   PA (4) 
and  PB(4)  have  a  target  in  common,  and  are  therefore  not 
feasible  for  the  two-aircraft  problem.   Since  no  other  ter- 
minal node  has  an  upper  bound  as  high  as  node  four,  the 
algorithm  branches  from  node  four  to  form  nodes  six  and 
seven. 

The  algorithm  continues  in  this  manner  with  the  upper 
bound  computed  for  each  new  node,  and  following  each  double 
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branching,  a  comparison  of  terminal  node  upper  bounds  and 
a  check  for  feasibility  is  performed.   When  nodes  eight  and 
nine  branch  from  node  three  and  the  bound  of  each  is  com- 
puted/ it  is  found  that  node  eight  has  the  highest  upper 
bound  of  terminal  nodes  five,  six,  seven,  eight,  and  nine. 
It  is  further  found  that  PA (8)  and  PB(8)  have  no  targets 
in  common.   Therefore,  the  algorithm  stops  and  PA (8)  and 
PB(8)  form  the  optimal  solution  with  a  value  of  26. 

Table  II  summarizes  the  progress  of  the  algorithm  in  the 
example.   Figure  12  illustrates  the  complete  branching  tree 
for  the  problem,  and  Figure  13  shows  the  optimal  path  of 
each  aircraft  through  the  target  area. 
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R(4)=A3,A6 


R(2)=A3 


R(3)=B- 


(9)=B1,B3 


♦Optimal 


Figure  12.   Branching  Tree  for  Test  Problem 
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PB(8) 


PA  (8) 


Figure  13.   Optimal  Path  of  Each  Aircraft 
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V.   VARIATIONS 

The  model  which  has  been  formulated  can  be  modified  to 
solve  more  difficult  problems  having  additional  constraints. 

One  constraint  would  be  to  limit  the  number  of  bombs  on 
each  aircraft,  thus  limiting  the  number  of  targets  allowed 
in  the  optimal  paths  of  the  aircraft.   To  incorporate  this 
restriction,  the  single  aircraft  D.P.  method  of  Section  I 
must  be  modified.   This  can  be  done  by  increasing  the  num- 
ber of  state  variables  from  one  to  two.   At  each  stage  there 
will  be  one  state  variable,  x  ,  representing  the  lateral 
position  of  the  aircraft  at  stage  n.   Another  state  variable, 
NBOMBS  ,  denotes  the  number  of  bombs  remaining  in  the  air- 
craft at  stage  n.   Although  the  computation  required  for  the 
D.P.  portion  of  the  algorithm  increases  exponentially  with 
the  number  of  state  variables,  practical  problems  can  still 
be  quickly  solved  with  this  added  constraint.   The  bomb 
limitation  has  been  incorporated  into  the  computer  program 
contained  in  Appendix  B. 

In  a  similar  manner,  a  restriction  on  the  total  lateral 
deviation  of  the  aircraft  or  on  the  total  number  of  course 
changes  allowed  in  the  target  area  could  be  considered. 

Another  modification  would  be  to  have  M  aircraft  attack- 
ing from  each  direction.   The  single  aircraft  D.P.  portion 
of  the  algorithm  becomes  instead  an  optimization  for  a  single 
group   of  aircraft,  where  a  separate  mutually  exclusive  path 
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is  computed  for  each  aircraft  in  the  group  (the  group  con- 
sisting of  all  aircraft  coming  from  one  direction) .   The  op- 
timization for  a  single  group  of  M  aircraft  is  a  dynamic 
programming  problem  with  M  state  variables.   The  lateral 
position  of  the  i   aircraft  at  stage  n  is  denoted  by  x.  . 
The  one  directional  problem  is  discussed  in  detail  in  Ref.  2. 
To  solve  the  two  directional  problem,  the  group  of  aircraft 
coming  from  one  direction  is  viewed  as  group  A,  the  other  as 
group  B.   Since  there  will  be  no  intersection  of  targets 
within  a  group,  the  only  concern  will  be  with  targets  in  com- 
mon between  the  two  groups.   As  in  Section  III,  the  highest 
valued  common  target  at  node  I,  XAB(I),  can  be  found  and  the 
branching  performed  with  one  node  of  the  branch  restricting 
group  A  to  exclude  target  xar(DA'  the  other  node  restricting 
group  B  to  exclude  XAB(DB- 

Again,  it  is  critical  to  realize  that  doubling  the  num- 
ber of  state  variables  far  more  than  doubles  the  computations 
required,  and  eventually  the  problem  will  become  computationally 
inf easible. 

The  next  modification  to  be  considered  is  the  problem 
of  aircraft  attacking  from  more  than  two  directions.   The 
general  theory  of  the  two-directional  problem  can  be  extended 
to  the  N  directional  case,  but  the  rules  governing  branching 
become  more  involved,  and  the  number  of  nodes  required  for 
solution  greatly  increases.   One  possible  approach  where  N  =  3 
will  be  briefly  considered,  with  the  aircraft  designated  A, 
B,  and  C. 

35 


At  each  node,  there  are  two  types  of  intersection  possible, 
a  two-aircraft  intersection,  and  a  three-aircraft  intersection. 
Branching  will  be  done  on  the  highest  valued  target  of  inter- 
section, whether  it  is  common  to  two  aircraft  or  three.   At 
node  I  this  point  will  be  designated  INT (I) .   Assume  that  at 
node  L,  INT(L)  =  A .  =  B . .   A  double  branch  would  emanate 
from  node  L,  one  restricting  aircraft  A  to  exclude  A.,  the  other 
restricting  aircraft  B  to  exclude  B..   Should  INT(L)  instead 

equal  A.  =  B.  =  C,  ,  a  triple  branching  would  be  required. 
1     3     K 

One  branch  would  restrict  aircraft  A  and  B  to  exclude  A.  and 

1 

B . ,  respectively.   The  second  branch  would  restrict  aircraft 
A  and  C  to  exclude  A.  and  C,  ,  respectively.   The  last  branch 
would  restrict  aircraft  B  and  C  to  exclude  B.  and  C,  ,  respec- 
tively. 

The  paths  at  each  node  are  calculated  as  described  in 
Section  I,  using  the  single  aircraft  D.P.  optimization,  sub- 
ject to  the  restrictions  above.   The  upper  bound  for  any  node 
is  the  summation  of  the  values  of  the  three  paths.   The  stop- 
ping condition  is  reached,  as  in  Section  III,  when  the  terminal 
node  with  the  highest  upper  bound  has  no  targets  in  common  on 
the  paths  of  the  three  aircraft. 

If  computer  time  is  critical,  a  suboptimal  solution,  as 
close  to  optimal  as  the  decision  maker  desires,  could  be  found. 
This  is  done  by  selecting  a  value  which  represents  the  maximum 
difference  the  decision  maker  can  tolerate  between  the  highest 
upper  bound  and  its  corresponding  lower  bound.   When  this  value 
is  achieved,  the  suboptimal  solution  is  obtained  by  randomly 
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assigning  the  targets  of  intersection  at  that  node  to  either 
aircraft,  thus  making  the  solution  feasible. 
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VI.   THE  COMPUTER  PROGRAM 

Appendix  B  contains  a  computer  program,  written  in  For- 
tran IV,  which  will  solve  the  target  optimization  problem  for 
two  aircraft  traversing,  from  different  directions,  a  target 
area  of  up  to  100  targets.   The  difference,  a,  in  the  initial 
courses  of  the  two  aircraft,  may  vary  between  0  and  3  60  degrees. 
The  aircraft  may  have  between  two  and  twenty  bombs  on  board. 
The  computer  program  gives  the  user  the  number  of  branchings 
required  for  solution,  the  optimal  value  of  the  targets 
chosen,  the  path  of  targets  each  aircraft  is  to  attack,  and 
a  plot  of  the  target  area  and  the  optimal  paths  through  it. 
Appendix  C  contains  the  output  from  a  one  hundred  target 
area,  with  an  allowable  course  deviation  of  forty-five  de- 
grees, six  bombs  per  aircraft,  and  a  =  ninety  degrees. 

The  input  parameters  are  the  total  number  of  targets, 
the  number  of  bombs  on  board  each  aircraft,  the  allowable 
course  deviation  angle,  the  difference  in  the  initial  courses 
of  the  aircraft,  and  the  location  and  value  of  each  target. 
The  target  positions  are  given  in  cartesian  coordinates,  with 
the  X  axis  perpendicular  to  the  initial  course  of  aircraft  A. 
The  angle  a  is  measured  counter-clockwise  from  aircraft  A, 
and  is  input  in  degrees.   The  allowable  course  deviation 
angle  is  identical  for  both  aircraft,  although  the  program 
could  easily  be  modified  to  allow  each  a  separate  deviation. 
This  angle  is  also  input  in  degrees. 
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The  program  begins  by  sorting  and  numbering  the  targets, 
first  with  respect  to  the  initial  course  of  aircraft  A,  and 
then  with  respect  to  aircraft  B. 

An  MN  by  MN  matrix  is  formed  for  determining  for  either 
aircraft  whether  one  target  position  may  be  feasibly  reached 
from  another.   A  "1"  indicates  feasibility,  and  a  "0"  infeasi- 
bility.   The  diagonal  elements  from  upper  right  to  lower  left 
are  all  zero,  indicating  that  the  aircraft  may  not  remain  at 
one  target  for  more  than  one  stage.   Denoting  any  element 
as  FEAS(I,J),  the  elements  below  the  diagonal  give  the  feasi- 
bility of  aircraft  A  going  from  target  I  to  target  J.   The 
elements  above  the  diagonal  give  the  feasibility  of  aircraft 
B  going  from  target  J  to  target  I.   This  matrix  eliminates 
the  need  to  geometrically  compute  feasibility  at  every  stage 
of  the  algorithm. 

A  vector  AB  is  formed  to  correlate  the  target  numbers  with 
respect  to  A  with  the  target  numbers  with  respect  to  B .   If 
AB(i)  ■  j  ,  then  A .  =  B . . 

The  path  restrictions  for  every  node  are  stored  in  a  100 
by  50  matrix  R.   The  matrix  permits  up  to  one  hundred  branch- 
ings of  the  algorithm,  and  restrictions  of  up  to  fifty  targets 
at  each  node.   Either  of  these  may  be  increased  by  the  user. 
The  restrictions  for  node  I  are  stored  in  row  I  of  R  by  denot- 
ing a  target  A.  as  negative  i,  and  a  target  B.  as  j,  thus 
signifying  whether  a  particular  numerical  element  restricts 
aircraft  A  or  aircraft  B. 
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The  single  aircraft  D.P.  portion  of  the  algorithm  is  per- 
formed in  subroutine  D.p.  of  the  program.  A  slightly  modified 
version  of  the  double  DO  loop  method  suggested  in  Ref .  2  is 
used.   It  is  presented  in  Fortran  in  a  simplified  form  in 
Figure  14. 

DO  30  I  =  2,MN 
DO  20  J  =  1,1-1 
C    Is  it  feasible  to  go  from  target  I  to  target  J? 

DO  10  K  =  2,NBOMBS 
C    If  I  have  K  bombs  on  board  at  target  I,  is  J 
the  best  target  to  go  to? 
10    Continue 
20    Continue 
3  0    Continue 
C    Trace  back  to  find  the  best  path 

Figure  14.   The  Simplified  Triple  DO  Loop  of 
Subroutine  D.P. 

A  one-bomb  limitation  is  not  allowed,  since  the  solution 
to  the  one-bomb  problem  is  merely  to  choose  the  two  highest 
valued  targets  and  assign  one  of  them  to  aircraft  A  and  one 
to  aircraft  B.   If  an  unlimited  number  of  targets  is  possible 
on  a  path,  as  might  be  the  case  in  planning  a  photo  recon- 
naissance mission,  the  bomb  limitation  should  be  completely 
removed  from  the  program,  rather  than  using  a  very  large 
number  for  the  limitation.   This  is  so  because  increasing  the 
number  of  bombs  on  board  increases  the  computational  effort 

required. 
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The  score  ordering  method  suggested  in  Ref .  8  was  tested 
in  the  program.   This  method  requires  that  for  target  I,  the 
targets  1-1,  ...,  1  are  stored  in  a  list  in  order  of  non- 
increasing  cumulated  value.   At  target  I,  the  list  is  scanned, 
starting  from  the  top,  until  a  feasible  target  J  is  found. 
The  cumulative  value  of  J  is  then  added  to  the  value  of  I, 
and  target  I  is  placed  in  the  list,  its  position  depending  on 
its  now  cumulative  value.   This  eliminates  the  need  to  scan 
all  the  lower  numbered  targets  from  target  I  to  find  the 
best  one.   It  did  in  fact  result  in  reduced  computation  for  up 
to  three  bombs  on  board.   But  with  more  than  three,  the  com- 
putation required  to  update  the  list  at  each  stage  outweighed 
the  savings,  and  therefore  score  ordering  was  not  included  in 
the  program. 

For  each  node  I  in  the  problem,  the  values  of  R(I) ,  UB(I), 
XAB(I),  and  VINT(I)  are  saved.   PA(I)  and  PB(I)  are  discarded 
as  soon  as  the  above  four  values  are  computed.   Since  the 
branching  tree  grows  horizontally  as  well  as  vertically  (see 
Figures  12  and  2  0) ,  a  large  amount  of  storage  would  be  used 
up  in  saving  the  paths.   For  this  reason,  new  paths  are  com- 
puted at  each  new  node  from  R(I)  and  X,B(I)  of  its  predecessor 
node  I . 

Figures  15,  16,  and  17  are  time  comparisons  for  various 
input  values  of  NBOMBS,  course  deviation  angles,  and  a.   The 
data  points  on  the  graphs  represent  averages  for  three  dif- 
ferent random  target  areas.   However,  the  trends  were  almost 
identical  for  each  set  of  targets. 
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Increasing  the  number  of  bombs  per  aircraft  causes  an  al- 
most linear  increase  in  execution  time.   An  increase  in  the 
course  deviation  angle  approximates  an  exponential  increase 
in  the  time  required.   Decreasing  a  causes  execution  time  to 
increase.   This  should  be  expected,  since  the  closer  to  parallel 
the  two  aircraft  are,  the  larger  the  number  of  targets  in  com- 
mon at  each  node  is  likely  to  be,  resulting  in  more  branching 
being  required. 

Figure  18  plots  execution  time  versus  the  number  of  tar- 
gets in  an  area.   Figure  19  illustrates  the  corresponding  num- 
ber of  branchings  required.   As  the  total  number  of  targets 
increases,  the  execution  times  tend  to  cycle.   One  possible 
explanation  for  this  is  that  as  the  number  of  targets  increases, 
the  amount  of  computation  in  the  dynamic  programming  subroutine 
increases.   However,  when  a  path  intersection  occurs,  more 
targets  provide  more  alternative  paths  which  may  be  feasible 
and  optimal,  thereby  reducing  the  number  of  branchings  re- 
quired.  For  certain  numbers  of  targets,  the  time  savings  from 
the  reduced  branchings  override  the  increased  subroutine 
computation  required,  thus  reducing  total  solution  time. 
Similar  trends  were  found  with  other  target  arrangements. 
This  again  points  out  the  unpredictability  of  the  computation 
required  for  solution. 
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Figure  16.   Execution  Time  Versus  Course  Deviation  Angle 
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Figure  18 .   Execution  Time  Versus  Number  of  Targets 
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VII.   CONCLUSIONS 

The  algorithm  and  computer  program  presented  can  solve  the 
two-aircraft  target  optimization  problem  for  up  to  one  hundred 
targets,  with  an  aircraft  payload  of  up  to  twenty  bombs,  using 
a  minimal  amount  of  computer  time.   In  addition,  it  can  be 
expanded  by  the  user  to  suit  his  specific  needs,  including 
more  targets,  larger  payloads,  more  aircraft,  and  a  multiple 
direction  attack  scenario. 

The  algorithm  itself  could  be  improved  if  there  were  some 
way  to  recognize  the  optimal  solution  before  the  stopping 
condition  was  met.   Oftentimes,  an  early  node  will  produce  the 
optimal  path,  but  the  algorithm  continues,  because  the  node's 
upper  bound  is  not  the  highest.   This  is  the  case  where  multi- 
ple optimal  solutions  exist,  and  each  level  of  the  branching 
tree  produces  equivalent  upper  bounds.   This  is  illustrated 
in  Figure  20.   The  upper  bounds  are  indicated  above  each  node. 
At  node  13 ,  an  upper  bound  of  thirty-seven  is  achieved  with 
no  targets  in  common  for  PA (13)  and  PB(13).   Yet  the  branch- 
ing must  continue,  since  other  terminal  nodes  have  higher 
upper  bounds.   Finally,  by  node  25,  it  is  realized  that  the 
solution  found  at  node  13  was  in  fact  optimal.   Had  this  been 
realized  at  node  13,  the  computation  required  could  have  been 
cut  in  half. 

The  computer  program  was  developed  to  test  the  algorithm, 
and  should  not  be  considered  as  an  end-product  software  package. 
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'Optimal 


Figure  20.   A  Branching  Tree  with  Alternative  Optimal 

Solutions. 


It  is  storage  inefficient,  especially  in  the  area  of  the  R  ma- 
trix and  the  fact  that  information  is  kept  in  storage  for  all 
nodes  rather  than  just  the  terminal  ones.   Improvements  in  these 
areas  could  be  implemented  if  desired. 

49 


APPENDIX  A 


GLOSSARY  OF  TERMS 


Allowable  Course  Deviation  =  the  allowable  number  of  degrees 

that  an  aircraft  is  permitted  to 
deviate  from  its  initial  course 
upon  entering  the  target  area. 

MN  =  the  total  number  of  targets  in  the  target  area. 

A.  =  the  i   target,  ordered  non- increasingly  by  its  parallel 

distance  from  the  boundary  at  which  aircraft  A  enters 
the  target  area,  i  =  1,  ...,  MN. 

B.  =  the   j        target,    ordered   non-increasingly   by   its  parallel 

distance  from  the  boundary  at  which  aircraft  B  enters 
the  target  area,  j  =  1,  ...,  MN. 
V[A.]  =  the  value  of  target  A- . 
V[B.]  =  the  value  of  target  B . . 

PA (I)  =  the  set  or  path  of  targets  for  aircraft  A  computed  at 
node  I.   It  is  feasible  and  optimal  for  the  single 
aircraft  problem. 
PB(I)  =  the  set  or  path  of  targets  for  aircraft  B  computed 

at  node  I.   It  is  feasible  and  optimal  for  the  single 
aircraft  problem. 
V[PA(I)]  =  the  summation  of  the  values  of  the  targets  in  PA(I) 
V[PB(I)]  =  the  summation  of  the  values  of  the  targets  in  PB(I) 
X^d)    =  the  target  common  to  both  PA  (I)  and  PB(I)  which 
has  the  highest  value. 
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X  (I).  =  XA_(I)  in  terms  of  aircraft  A. 
Xjgdlp  =  X^d)  in  terms  of  aircraft  B. 
VINT (I)  =  the  summation  of  the  values  of  all  targets  common 

to  both  PA(I)  and  PB(I) . 
UB(I)    =  V[PA(I)]  +  V[PB(I)]  =  an  upper  bound  on  the  optimal 

solution  at  node  I. 
LB (I)    =  UB(I)  -  VINT (I)  =  a  lower  bound  on  the  optimal 

solution  at  node  I . 
R(I)     =  a  set  of  restrictions  on  PA(I)  and  PB(I)  consisting 

of  a  listing  of  targets  which  each  must  exclude. 
NBOMBS   =  the  number  of  bombs  with  which  an  aircraft  enters 

the  target  area. 
a  =  the  difference  in  the  initial  courses  of  the  two  aircraft, 
measured  in  degrees  counter-clockwise  from  aircraft  A. 
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APPENDIX    B 


THE    COMPUTER   CODE 


C    A   TWO    DIRECTIONAL    TARGET    OPTIMIZATION    MCDEL    USING    THE 
C    BRANCH    AND    BOUND    ALGORITHM    WITH    DYNAMIC    PROGRAMMING. 
C 

c 

C    THIS     PRCGRAM    COMPUTES    OPTIMAL    TARGET    PATHS    FOR    TWH 

C    AIRCRAFT    TRAVERSING    A    TARGET    AREA    FROM    CIFFEREMT 

C    DIRECTIONS.     IT    WILL    PERFORM    UP    TO    100    BRANCHINGS.    THE 

C    TARGET    POSITIONS    ARE    IN    TERMS    OF    THE    CARTESIAN    PLANE,    WITH 

C    THE    X    AXIS    PERPENDICULAR    TO    ONE    OF    THE    AIRCRAFT'S     INITIAL 

C    COURSES.    THE    NUMBER    OF     eOMBS     UN    BOARD    EACH    AIRCRAFT    MUST 

C    BE    BETWEEN    TWO    AND   TWENTY.    THE   TARGET     AREA    MAY     HAVE    UP    TO 

C    ONE    HUNCRED    TARGETS. 

C 

C    THE     INPUT    PARAMETERS    ARE    : 

C      MN    =   TOTAL    NUMBER    OF    TARGETS 

C       NB0V3S    =    NUMBER    OF    BOMBS    ON    BOARD    EACH    AIRCRAFT 

C      ANGLE    =    ALLOWABLE    COURSE    OEVIATION    (     IN    DEGREES     ) 

C      ALPHA    =    DIFFERENCE    IN     INITIAL    COURSES    OF    THE    TWO    AIRCRAFT 

C       AX(I)     =    X    COCRDINATE    OF    TARGET     I 

C      AY(I)     =    Y    COORDINATE    OF    TARGET    I 

C      VA(  I)     =    VALUE    OF    TARGET     I 

C 

C    MAIN    PROGRAM 

C 

C 

C    DIMENSION    THE    ARRAYS    AND    ASSIGN    VARIABLE    TYPE 

REAL    AX  (100),  AY  (100  ) , V  A ( 100 ) , XA (  10Q),YA(100)  ,  BX(IOO)  , 
1BYU00)  ,VBUO0)  ,XB(100),YB(100),TVA(100  ),TVB(100), 
1 VI  NT (100) ,UB( 100) , V(100) 

INTEGER    FEAS<100,100  ),R(100,50)  ,AB(  100),XA8(  ICO), 
IPATH(IOO)  ,PA(  100) 
INTEGER    COUNT 
C 
C 

C    READ    IN    THE    NUMBER    OF    TARGETS,     THE    NUMBER    OF    BOMBS    PER 
C    AIRCRAFT,    THE    DEVIATION    ANGLE,     AND   THE    CIFFERENCE    IN 
C    INITIAL    COURSES    OF    THE    TWO    AIRCRAFT. 
READ     (5,10)     MN,  NBOMBS,  ANGLE, ALPHA 
10    FORMAT    (2I4f2F5.0) 
C 
C 

C    READ    IN    THE    TARGET   POSITIONS    AND    VALUES 
00    20    1=1, MN 

REAC    (5,30)     AX(I)  ,  AY  ( I  ),VA(I) 
20    CONTINUE 
30     FORMAT     (3F12.0) 
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c 
c 

C    SORT    TARGETS    FOR    AIRCRAFT    A 

MNM1    =    MN-1 

DO    5  0    I»i,MNM-l 

IP1    =    1+1 

IOPT   =    0 

CO    40    J=IP1,MN 

IF     (AY(I)  .GE.AY(  J)  )    GO    TO   40 

TEMPY    =    AY(  I  ) 

TEMFX    =    AX  (I) 

TEMPV    =   VAU  ) 

AY(I)    =    AY(J  ) 

AX(  I)    =    AX(J1 

VA(I)     =    VA(J) 

AY  (J)    =    TEMPY 

AX(J)     =   TEMPX 

VA(J)    =    TEMPV 
40    CONTINUE 
50    CONTINUE 
C 
C 
C    SORT    TARGETS    FOR    AIRCRAFT    B 

ALPHA    =    ALPHA*3. 141592654/180.0 

DO    60    I=L  ,MN 

3X(I)     =    AX(  I)  *COS(  ALPHA)+AY(I  )*SIN(ALPHA) 

8Y(I)    =    -AX  (  I  )*S  IN(  ALPHA)+AY(  I)*COS(ALPHA  ) 

V8(I)     =    VAU  ) 
60    CONTINUE 

CO    90    1=1  ,MNM1 

IP1    =    IH 

DO    70    J=IP1,MN 

IF     (BY(I)  ,GE.3Y( J) )    GO   TO    70 

TEMPY    =   3  Yd  ) 

TEMPX    =    BX(I  ) 

TEMPV    =   V8(I  ) 

BY(I)    =    BY( J) 

BX(I  )    =    BX(J  ) 

V8( II     =    VB{ J) 

3Y(J)    =    TEMPY 

eX  (J)    =    TEMPX 

VB(J)    =    TEMPV 
70    CONTINUE 
80    CONTINUE 
C 
C 
C    COMPUTE    FEASIBILITY    MATRIX 

ANGLE    =    (  (90.  0-ANGLE)*3. 141592654/180. 0)-1.0E-8 

FEAS  (1,1  )    =    0 

00    120    1=2, MN 

IM1    =    1-1 

FEAS(Irl)    =    0 

DO    110    J=1,IM1 

IF    (  AX(  I  )  .EQ.AX(  J  )  )    GO    TO    90 

IF     (BXU)  .EQ.BXUH     GO    TO    100 

SLOPEA    =    ATAN(ABS( ( AY ( J)-AY(I ) ) /( AX( J)-AX(  I)  )  )  ) 

SLOPEB    =     AT  AN  (ABS  (  (  BY(  J  )-BY(  I  )  )/(BX(  J)-BX(  I)  )  )  ) 

FEASd  ,J)     =    0 

PEAS  (J,  I)     =    0 

IF     (SLOPEA. GT.AMGLE)    FEAS(I,J)  =  1 

IF    (SLOPEB. GT  .ANGLE)    FEAS(J,I)=1 

GO    TO    110 
90    FEASd, J)    =    1 

GO    TO    110 
100    FEAS(J,I)    =    1 
110    CONTINUE 
120    CONTINUE 
C 

c 


" 


C    SET    UP    AB    CORRELATION    VECTOR 
DO    140    L=1,MN 
DO    130    M=1,MN 
TWOPIE    =    2.0*3.141592654 
TPPA    =   TWCPIE-ALPHA 

IP     (ABS( AX(L)-(BX(M)*COS(TPMA)+BY(M)*SIN(TPMA)  ) ) . GT . 
11.0E-2)    GO    TO    130 

IF     (ABS(  AY(L)-(-BX(  V  )*S  IM(T  PM  A)  +BY  (  M  )*COS  (  T  PM£ )  )  )  .GT 
11.0E-2)    GO    TO    130 
48  (L )    =    M 
GO    TO    140 
130    CONTINUE 
140    CONTINUE 
C 
C 

C    BEGIN    BRANCH    ANO    BOUND 
C 
C 

C    ZERO    OUT    RESTRICTION    MATRIX 
DO    160   L=l, IOC 
CO    150    M=l,50 
R(L,M)     =    0 
150    CONTINUE 
160    CONTINUE 
C 
C 

C    DETERMINE    R    VECTOR    FOR    NODE    I 
XABd  )    =    0 
I    =    0 
J    =    1 
K    =    0 
170    I    =    1+1 

IF     (I.E0.101)     GO    TO    490 
DO    130    L=l,50 
LAST    =    L 

IF    (RU,L  ).EQ.O)    GO    TO    190 
IF     (L.EQ.50)    GO   TO    510 
R(I  ,L)    =    R(J,L) 
180    CONTINUE 
190    IF     (I.EQ.21     K    =    0 

IF    (K.5Q.1)    GO    TO    200 
R(I,LAST)     =    -XAB(  J  ) 
K    —    1 
GO    TO    210 
200    R(I  ,LAST)     =    AE(X  AB(  J)  ) 
K    =    Q 
C 
C 

C    ADJUST    TARGET    VALUES    TO    EXCLUDE    TARGETS    IN    RU) 
210    DO    230    L=l,  50 

IF    (R( I,L) .EQ.O)    GO    TO    240 
IF    (R(I,L  ).GT.O)    GO    TO    220 
TVA(-R( I,LJ  )    =   VA(-R(I,L) > 
VA(-R(I ,L) )     =   -10000.0 
GO    TO   23  0 
220    TV8(R(I,D)     =    VB(R(I»L)J 

V3(  R{  I,L)  )    =    -10000.0 
230    CONTINUE 
C 
C 

C    CALL    DP    FOR    A    AND    B 
240    ICHK    =    1 

CALL    DP     (FEAS,VA,ICHK,MN,NBOMBS, PATH, W, COUNT) 
VPA    =    W 
NPA    =    COUNT 
00    250    L=1,NPA 
PA(L)     =    PATH(L) 
250    CONTINUE 
IChK    =    2 

CALL    DP    (  FEAS,VB,ICHK,MN,NBOMBS,  PATh,W, COUNT  ) 
IF    (  IOPT.EO.l  )    GO    TO    350 
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c 

C    RESET    TARGET    VALUES 

DO    270    L=l,50 

IF     (R(I,L).EQ.O)    GO    TO    280 

IF     <RU,L  I.GT-.O)    GO    TO    260 

VA(-R(I»D)     =    TVA(-R(I  ,L)  ) 

GO    TO    270 
260    VBIRdtU  )    =    TVB(R(I,L)  ) 
270    CONTINUE 
C 
C 

C    COMPUTE    XAB(I),    VINT(I),    AND    UBU) 
280    XMAX    =    0.0 

SUM    =    0.0 

DO    300   L=1,NPA 

CO    290    M=l,  COUNT 

IF     (AB(PAU)  J.NE.PATH(M)  )    GO    TO    290 

SUM    =    SUM+VA(PA(L) ) 

IF     (VA(PA(L)l  .LE.XMAX)     GO    TO    300 

XMAX    =    VA(PA(D) 

XAE(I)    =    PA(L  ) 

GO    TC   300 
290    CONTINUE 
300    CONTINUE 

VINT!  II     =    SUM 

UBU)    =    VPA+W 

JF     U.EQ.l.AND.VINTd  )  .  EO  .0  .0  )    GO    TO    340 

IF     (K.EQ.  1)     GO    TO    170 

UBU  )    =    0.0 

c 
c 

C    CHOOSE    FIGHEST    UPPER    BOUND    OF    ALL    TERMINAL    NODES 
C    AND    CHECK    IF    IT    IS    FEASIBLE 

UBMAX  =  0  .0 

DO  310  L=lf  I. 

IF  (UB( L)  .LE. UBMAX)  GO  TO  310 

UBMAX  =  UB(L  ) 

J    =    L 
310    CONTINUE 

IF     (VINT  (J  J.EG.O.O)     GO    TO    330 
C 
C 

C    CHECK    FCR   TIE    FOR    HIGHEST    UPPER    BOUND. 
C    IF    THERE    IS    A    TIE,    IS    IT    FEASIBLE? 

DO    320    L  =  lt  I 

IF     (UBU)  .NE.UBU  )  )    GO    TO    320 

IF    (  VINTl  D.NE.O.O)     GO    TO    320 

J    =    L 

GO  TC  330 
320  CONTINUE 

GG  TO  170 
C 
C 

C    IF    NECESSARY,     RECOMPUTE    OPTIMAL    PATH 
330    NODES    =    I 

IF    (  I.EQ.J  )    GC   TO    350 

I    =    J 

IOPT    =   1 

GO    TO    210 
C 
C 
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C    PRINT    OUT    THE    NUMBER    OF     BRANCHING    NODES,    THE    OPTIMAL 
C    VALUE,    AND    THE    OPTIMAL    PATHS. 
340    NODES    =    I 
350    WRITE    (6,360) 
36C    FORMAT    CI') 

WRITE    (6,370)    NODES 
370    FORNAT    ('     OPTIMAL    SOLUTION    FOUND    IN', 13, •     BRANCHINGS', 
1////////) 
WRITE     (6,  380)    UB(  J  ) 
380    FORMAT    ('     THE   VALUE    OF    THE    OPTIMAL    SOLUTION     IS     »,F12.4 
It'////////) 
WRITE    (6,390) 
390    FaR^AT    ('     OPTIMAL    AIRCRAFT    PATHS',///) 

WRITE    (6,400) 
400    FORMAT     (6X.' AIRCRAFT      A',///) 

WRITE    (6,410) 
410     FORMAT    (  4X,  •  X  •  ,12X,  •  Y  • , //) 
DO    420    1=1, NPA 
XA(I)    =    AX(  PA(I)  ) 
YA(I)    =     AY  (  P  A  (  I )  ) 
WRITE    (6,430)    XAU  ),YA(  I) 
420    CONTINUE 
430    FORMAT     (IX,  F8  .2,  5X,  F8  .2,/ ) 

WRITE    ( 6,440) 
440    FORMAT    (///,  6X ,' AIRCRAFT      B',///J 
WRITE    (6,410) 
00    450    1  =  1, COUNT 

XB(I)    =     BX(PATH(  I)  )  *COS  (TPMA)+6Y(  PATH(  I  )  )*SIN  (TPMA) 
YB(I)     =    -BX(  PATH(  I  )  )*SIN(TPMA)  +  EY<P4TH(  I)  )*COS(TPMA  ) 
WRITE    (6, 430)     XB( I)  ,YB( I) 
450     CONTINUE 
C 
C 

C    PLOT    THE    TARGETS    AND   THE    PATHS 
WRITE    (6,460) 
460    FORMAT     CI') 

CALL    PLOT  P    (AX, AY,  MN,1  ) 
CALL    PLOTP     (XA,YA,NPA,2 ) 
CALL    PLOTP     (Xe,YB,  COUNT  ,3  ) 
WRITE    ( 6,470) 
470    FORMAT    ( / //, 40X, ' P ATH    OF    AIRCRAFT    A    DENOTED    BY    +• ) 

WRITE    (6,480) 
48C    FORMAT    (/,4DX,'PATH    OF     AIRCRAFT    8    DENOTED    BY    *»  , 
1////////////) 
GO    TO    530 
C 
C 

C    CHECK    FOR    DIMENSIONAL    ERRORS 
490     WR  ITE    (  6,500) 

500    FORMAT    (////,'    REQUIRED    BRANCHING    EXCEEDS    100    NODES', 
1////) 
GO    TO    530 
510    WRITE    (6,520) 

520    FORMAT    (////,«    RESTRICTION    VECTOR    EXCEEDS    50',////) 
530    STOP 
END 
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SUBROUTINE    DP     (FEAS, V,  ICHK,MN ,NBOMBS , PATH, W, C CUNT) 
C 

C    THIS    SUBROUTINE    COMPUTES    THE    OPTIMAL    PATH    FOR    A    SINGLE 
C    AIRCRAFT    TRAVERSING   THE    TARGET    AREA. 
C 

c 

REAL    VV(I00,20),  V(  100)  ,CVAL(20> 

INTEGER    P  (100  ,20),  FEAS  (100, 100),  MOVE  (20),  PATH  (100) 
INTEGER    CCUNT 
CVAL(l)    =    0..0 
NOVEU)     =    200 
DO     20    K=1,NB0MBS 
VV  <1,K)    =    V(l) 
P(  1,K)    =    200 
20    CONTINUE 
C 

c 

C    I     IS    THE    TARGET    THAT    THE   AIRCRAFT    IS    AT    PRESENTLY 

DO   7  0    1=2  ,MN 

W(  1,1)    =    V(  I  ) 

CO    30    K=2,NB0MBS 

CVAL(K)     =    0.0 

MOVE(K)    =    200 

VV  (If Kl    =    V(  I  ) 
30    CONTINUE 

IM1    =    1-1 
C 
C 
C    J    IS    THE    TARGET    THAT    THE    AIRCRAFT    MAY    GO    to    FROM    I 

00    50    J=1,IM! 

IF    (  ICHK.EQ.l.AND.  FEAS(  I,  J)  .EQ.O)    GC    TO    50 

IF     (ICHK.EQ.2  .ANO.FEAS(  J,  I  )  .EQ.O)    GOTO    50 
C 
C 
C    K    IS    THE    NUMBER    OF    BOMBS    LEFT    AT    TARGET    I 

DO    40    K=2,NB0MBS 

IF     (( VV( I,K  )  +  VV(J,(K-l)  ) ).LE.CVAL(K)  )    GO    TO    4C 

CVAL(K)    =    VV(  I,K)+VV  (J  ,  (K-l)) 

MOVE(K)     =    J 
40    CONTINUE 
50    CONTINUE 

00     60    K=1,NB0MBS 

IF     (CVAL(K)  .NE.0.0  )    VV(  I  ,K  )  =C  VAL  (  K  ) 

P(I,K)     =    MCVE(K) 
60    CONTINUE 
70    CONTINUE 
C 
C 
C    COMPUTE    FIRST    TARGET 

CVAL(NBOMBS)     =    0.0 

DO    80    J=1,MN 

IF     (VVUfNBCMBS  J.LE.CV  AL(N80MBS  )  )     GOTO    30 

CVAL(NBOMBS)     =    VV(J,NBOMBS) 

MOV  E(N  BOMBS  )     =    J 
80    CONTINUE 

W    =    CVAL(  NBOMBS) 
C 
C 
C    COMPUTE    OPTIMAL    PATH 

PATh(l)    =    MOVE(NBOMeS) 

COUNT    =    1 

LEFT    =    NBOMBS+1 

CO    90    I=2»MN 

COUNT    =    I 

LEFT    =    LEFT-1 

PATH(I)    =    P(  (  PATH(I-l)  ),LEFT) 

IF     (P((PATH( I )), (LEFT-1 )).EQ. 200)     RETURN 
90    CONTINUE 

RETURN 

END 
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APPENDIX   C 


THE    COMPUTER    OUTPUT 


OPTIMAL    SGLUTION    FOUNO     IN       7    6PANCKINGS 


THE    VALIE    GF    THE    CPTIfAl    SCLLTICN     IS  216.2231 


OPTIMAL  AIRCRAFT  PATHS 
AIRCRAFT   A 


47.71  48.91 

31.71  S5.13 

171. SI  297.32 

317. C5  445.82 

6d9.03  886.22 

754.56  585.23 

AIRCRAFT   6 


794.25  411.59 

4S1.19  389.44 

424.49  262.81 

363.56  378.05 

328. 59  224.02 

209.62  363.97 
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